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MODAL VECTOR ESTIMATION FOR CLOSELY-SPACED- FREQUENCY MODES 


1 , Introduction 

In describing the dynamical behavior of a complex structure modal 
parameters are often employed; undamped natural frequency, mode shape, 
modal mass, modal stiffness, and modal damping. From both an analytical 
standpoint and an experimental standpoint, identification of modal para- 
meters is ma^'te more difficult if the system has repeated frequencies or 
even closely-spaced frequencies. The more complex the structure, the 
more liLely it is to nave closely-spaced frequencies. In many cases this 
fact makes it difficult to determine valid mode shapes using single-shaker 
test methods. By employing band selectable analysis (zoom) techniques 
and by employing Kennedy-Pancu circle fitting or some multiple degree of 
freedom (MDOF) curve-fit procedure, the usefulness of the single-shaker 
approach can be extended. However, for many structures such procedures 
may still not be sufficient to give accurate modal representations. 

It is the purpose of this paper to discuss techniques for obtaining 
improved modal vector estimates for systems with closely-spaced-frequency 
modes. 

2. Cl osely-Spaced-Frequency-Modes 

It is well-known that systems with unique frequencies possess unique 
mode shsipes, while systems with repeated frequencies do not possess unique 
modes corresponding to the repeated frequencies, but rather, they possess 
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subspaces of modes. This Is Illustrated by the simple 2 degree-of-freedom 
(DOF) systems of Figure 1. Figure 2(a) shows the frequency response func- 
tions (R^F's) of a 2D0F system with w1 del y-separated-frequency modes 
(5,0Hz and 5,5Hz), while Figure 2(b) shows a similar system with closely- 
spaced-frequency modes (5.00Hz and 5.05Hz). 

The frequency spacing of the 2D0F system In Figure 2 Is controlled 

n ) 

by the strength of the coupling spring, Klosterman. * made a thorough 
investigvntlon of the dynamics of weakly coupled systems, such as the 
types illustrated In Figure 3, and concluded that, for such systems, 

"close agreement between computed and experimentally measured mode shapes 
Ci-/'‘not be expected." This results from the fact that "coupled system 
mode shapes are very sensitive to small variations In the subsystems." 

Table 1 illustrates the effect of small variations in the mass of one 
subsystem on the coupled system modes of a 2D0F spring-mass oscillator 
when the coupling is "weak." It can be seen that, as noted by Klosterman, 
small changes In the subsystem properties have an enormous effect on the 
system modes when the system consists of weakly coupled subsystems. 

Figures 4 and 5 show FRF's of a moderately complex piece of space 
hardware. Figure 4 shows the Bode plot of a drive-point FRF over the full 
data acquisition bandwidth, while Figures 5(a) and 5(b) show two Argand 
plots over a 10-30H2 bandwidth. Clearly, this structure has closely- 
spaced-frequency modes. In actual fact, the structure is comprised of 
a tubular framework supporting several virtually identical honeycomb panels. 
Thus, this structure is a classic example of a system which consists of 
weakly- coupled subsystems. 
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Although Klost'/irman has noted that getting experimentally determined 
modes to agree with analytically determined modes may be a hopeless under- 
taking, we will address the slightly different proposition of determining 
"good" (In some sense) experimental modes for systems having closely- 
spaced frequencies. 

3 . Modal Vector Estimation Using a Single Column/Row of the Frequency 

Response Function Matrix 

Techniques for estimating experimental mode shapes by using a single 
column or a single row of the FRF matrix are well known and are Incorporated 
In various modal analysis software packages such as MODAL-PLUS. Appendix A 
gives the mathematical basis for using a single row/column of the FRF matrix 
In identifying modal parameters. The most frequently employed techniques 
for estimating medal vectors from FRF's are the quadrature response method 
and the Kennedy-Pancu circle-fit method. Figure 6(a) shows the 
correct modes of the 2D0F systems of Figures 2(a) and 2(b) along with the 
quadrature response "experimental" modes which would be Inferred from fre- 
quency response functions. When the frequencies are not closely spaced, 
acceptable results are obtained by the quadrature response technique, as 
seen in Figure 6(b). However, when the frequencies are closely spaced, 
there is difficulty, first, In identifying the natural frequencies at which 
to Identify the modes; and modes based on one shaker location do not agree 
with modes based on another shaker location, as can be seen by comparing 
Figures 6(c) and 6(d) . 
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The use of quadrature response techniques in situations where fre- 
quencies are closely spaced clearly leads to erroneous modal vector esti- 
mates, It could be shown that even circle-fitting and MDOF curve fitting 
would fail to produce good modal vector estimates in cases such as this, 
Because of the difficulties in obtaining accurate modal vectors 
for sys'tems with closely-spaced -frequency modes by using only a single 
column or row of the FRF matrix, it becomes necessary to examine the possi- 
bility of using several rows/columns of the FRF matrix to obtain improved 
modal vector estimates. 

4. Modal Vector Estimation Using Several Columns/Rows of the Frequency 
Response Function Matrix 

Figure 6 illustrates the difficulty in obtaining valid mode shapes 
from single-shaker data when the system has dosely-spaced-frequency modes. 
If only one shaker location is employed, some modes may not be excited at 
all, and there may be superposition of contributions of several dosely- 
spaced-frequency modes leading to inaccurate modal estimates. If additional 
shaker locations are employed, there may be ambiguity as to which shaker 
location gives the best estimate of a certain mode. Several techniques 
have been employed to combine the information from several rows/columns 
of the FRF matrix. These include: 

Modal Vector Weighting 

Symmetry/Antisymmetry Weighting 
Analytical Mode Weighting^^^ 
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Ad Hoc Sort Logic Averaglng^^^ 

Modal Tuning; Force Appropriation 
Asher Method^®”^^ 

Extended Asher Method 

Minimum Coincident Response Method^^^^ 

Multishaker Sine Dwell Testlng^^^^ 

Figures 7(a) and 7(b) Illustrate the problems caused by the super- 
position of the effects of closely-spaced-frequency modes. In Figure 7(a) 
It 1s seen that the presence of two modes is masked, and It appears that 
a single mode with greater amplitude and greater damping, or possibly even 
nonlinear effects, is present. Figure 7(b) shows the effect of mode 
superposition on a non-driving-point FRF. Here two modes are clearly 
evident, but the apparent amplitude of each mode Is reduced from its true 
value, The methods listed above may be employed In an attempt to separate 
out the contributions of Individual modes In such situations. 

Figure 8 illustrates the basis of symmetry/antisymmetry averaging. 
Ground vibration testing of aircraft frequently employs symmetric/anti- 
symmetric shaker configurations because of the inherent synwetry of the 
structure. Even so, thf% • sy be many closely-spaced-frequency modes 
In either the symmetric or the antisymmetric set of modes of an aircraft. 

From Equations (A-7) and (A-10) It can be seen that, when the system 
has distinct roots, pj^ , every row and column of the residue matrix 
aj^ contains the modal vector U[^ multiplied by a component of itself. 
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Richardson and Kniskern^^^ extended the sort logic employed in producing 
mode shapes from a single row or column of the residue matrix and proposed 
a sort logic which would make use of data from multiple rows/columns of 
the residue matrix. Although they demonstrated some improvement in 
modal vector estimation, they recognized that their sort algorithm was 
based on “rather arbitrary rules," and that other "more 'optimum'" 
approaches might be possible* 

Modal tuning, or force appropriation, is a third approach to employ- 
ing multiple columns or rows of the FRF matrix to obtain system modes. 

The name "Asher method" is frequently applied to one such version of 
modal tuning. Extensive research has been done on this method in France, 
where the technique is referred to as "force appropriation." Reference 9 
is just one of a number of the French papers dealing with the use of 
force appropriation for modal vector estimation. 

The Asher method of modal tuning is based on the fact that when 
a system is excited at an undamped natural frequency by the appropriate 
force vector F^, » the response will be in quadrature with the excita- 
tion, and the quadrature response will correspond to the undamped natural 
mode <p^ That is, for monophase harmonic excitation of an nDOF 

system, 

X = H F ^ (Hr + iHj)F (1) 

nxl nxn nxl 

For the response to be in quadrature with the excitation 
HrF = 


0 


( 2 ) 
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Thus, the natural frequencies are determined from the condition that 

det(Hf^(( 0 j,)) • 0 , r * 1, 2, .... n (3) 

The force vector required to tune the mode Is given by Equation (2). 
Finally, the mode corresponding to the frequency Is given by 

X(o)j.) * 1Hj(o)^)F^ (4a) 

or 

Figure 9 shows the mechanics of how Asher tuning employs two columns of 
a FRF matrix to tune a moci^. 

in Equations (1) through (4) and In the example presented 1n 
Figure 9, It was assumed that all columns of the FRF matrix were available, 
1,e. that there could be forces Imposed at all points where displacement 
(velocity, acceleration) measurements were made. In actual fact, the 
number of possible excitation points, p , 1s usually one or two orders 
of magnitude smaller than the number of response points, n . In this 
case, Equations (1) through (4) are modified as follows: 

X = H F p<n (5) 

pxl pxp pxl 

F = 0 (6) 

tHere it will be assumed that n , the number of response points, 1s 
sufficiently large to constitute the "order of the system" as far as 
the behavior of the system in the frequency range of Interest Is concerned. 
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% 

A A 

det(HR({o)) * 0 

(7) 


A /V AAA A 

X (Up) • 1Hj(Up)Fp H- 0p 

(8) 


In Reference 7 It is shown that the frequencies and mode shapes obtained 
by applying Equations (5) through (8) to a pxp subset of the FRF matrix 
are frequently good approximations to the true frequencies and modes, 
but that “spurious modes” may also be produced w'fen p<n . Modal tuning 
according to Equations (5) through (8) will be referred to as “standard 
Asher" (SA) tuning. Examples of standard Asher tuning will be presented 
in Section 5. 

As noted, when the number of excitation points is small, i.e, 
p«n , the use of standard Asher tuning may lead to "spurious modes." 
These will have a response in quadrature with the excitation at the p 
excitation points, but at other measurement points the response may 
have a significant coincident component, thus giving the appearance 
of a "complex mode." Several approaches have been suggested for employ- 
ing a non-square submatrix of the FRF matrix. ^ The technique 

fill 

proposed by Ensrainger and Turner' ' may be referred to as the “minimum 
coincident response (MCR) method." Thus, for harmonic excitation, 

~ ~ ~ ~ p«n 

X = H F = (Hr + 1Hj)F (9) 

mxl mxp pxl m<n 


Xr = HrF 


(10) 
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Ut the coincident response norm be defined by 

e ■ x/Xr (11) 

Then, the force vector F Is appropriated such that c Is minimized 

0i0 

subject to an amplitude constraint. That Is, F(o)) and e(w) are 
determined as functions of o) only by determining an expression for 

A 

the F{(o) which minimizes e subject to an amplitude constraint. 


Min e (F(o))) ) 



m 1 


/ 


F(o)) , e(m) 


( 12 ) 


where X^j Is the quadrature response at the 1th response point. This 
results In a force vector given by 


F 


TiTPrr'T ^R^ ^11 




(13) 


where Hu Is the ith row of the Imaginary part of the FRF matrix. 
Finally, the frequencies are selected to minimize e(w) , Equation 
(13) Is used to compute the corresponding F^ , and the mode shape Is 
determined from the quadrature response 


Xj = 


(14) 


10 


"I 
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Ensmlnger and Turner reported good results from the application 
of the minimum coincident response method to synthesized froquency 
response functions based on ground vibration test data obtained for the 
mated 747-$pace Shuttle Orbiter, They also applied the standard Asher 
method, Equations (5) through (8), to the same synthesized FRF's. They 
concluded that "this approach (1.e, the minimum coincident response 
method) was found to provide more consistent results than that based 
on direct determination of roots of the real part of the transfer function 
matrix {i.e, the standard Asher method)." Gold and Hallauer^®^ discussed 
the application of the standard Asher method to synthesized '^RF's, but 


W*1 A 

ultaui c 


to successfully apply the method to actual experimental data 


5, Modal Tuning Applications 

A research program is being conducted to evaluate the use of the 
standard Asher method and the minimum coincident response method In 
estimating modal vectors of systems having closely-spaced-frequency 
modes. The specific questions being addressed are; 

(1) Is modal tuning using either the standard Asher(SA) method 
or the minimum coincident response(MCR) method effective in producing 
modal vector estimates in situations where the results of single shaker 
methods are not adequate? 

(2) How do the modal vector results obtained by the (SA) method 
and the (MCR) method compare? 

(3) What is the effect of the frequency resolution used in the 
synthesis of FRF's on the modes produced by the SA method and the MCR 
method? 


n 


% 




(4) How can residuals be employed In the synthesized FRF’s 
used as Input to the modal tuning methods? 


Although work is still In progress on this research project, some 
tentative results In response to questions (1) through (3) will be pre- 
sented. Figure 10 shows a flow chart of the major steps employed In the 
modal tuning studies. 

Figure 11 shows a coupled-beam model which permitted experimental 
data to be acquired for a system with closely-spaced-frequency modes. 
Figure 12 shows a typical FRF and the results of using the MDOF curve- 


fit procedure In MODAL-PLUS (the GE and GA commands). The original 


experimental bandwidth was 10 Hz, which concentrated the spectrum around 
the lowest pair of closely-spaced-frequency modes. Table 2 gives the 
MDOF curve-fit results for the eight FRF's (2 reference points, 4 response 
points). From Table 2 the following may be noted; 

(a) Each FRF has two dominant roots In the 118-121 Hz range. 

(b) Some of these dominant roots have phase angles significantly 
different than +90®, 

(c) The "modes" based on having a single shaker at 1 appear to be 
at approximately 119.0 Hz and 119.5 Hz; while the "modes" based on having 
a single shaker at 6 appear to be at 118.8 Hz and 120.0 Hz, However, the 
119.0 Hz "mode" Is radically different than the 118,8 Hz "mode," 


The MDOF results presented In Table 2 were used as input to a 
program which first synthesizes analytical FRF's and then employs the 
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synthesized FRF's as Input to the SA and MCR modal tuning algorithms. 

The modal tuning results for bandwidths of 10 Hz and 0.256 Hz are shown 
In Figures 13 and 14 respectively, and the results for several bandwidths 
are given In Table 3, From Table 3 the following may be noted; 

(a) Two modes at approximately 118,8 Hz and 119.0 Hz are consis- 
tently Identified by both methods. 

(b) Reducing the analysis bandwidth reduces the phase error of 
the identified modes. 

(c) The mode at 118,8 Hz appears to be much better than the mode 
at 119,0 Hz, I.e. there appears to be less phase error. However, it 
should be noted that the large phase errors In mode 2 are associated 
with relatively small amplitudes. In Figures 13(b) and H(b) it may be 
noted that the minimum of e(w) is much less sharp for the 119.0 Hz 
mode than for the 118.8 Hz mode. This sharpness may be related to the 
phase coherence of the mode, but this relationship has not been thoroughly 
eval uated. 

(d) For the 118.8 Hz mode there is good agreement between the SA 
mode and the MCR modi?. For the second mode the agreement Is not so good. 

(e) The tuned modes do not exhibit the symmetrical antisymmetric 
behavior which was expected because of the apparent symnetry of the 
structure. The explanation for this discrepancy may lie in the mode 
sensitivity of weakly coupled systems as previously discussed. 

Attempts to employ the SA modal tuning techr^..Me on some of the data 
from the structure whose FRF's are shown in Figures 4 and 5 indicated that 
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1t would be necessary to include residuals in the synthesized FRF's. Work 
on residuals is in progress and no results are yet available, 

6, Conclusions 

On the basis of the research described above, the following con- 
clusions may be stated: 

(a) Standard Asher tuning and minimum coincident response tuning 
both provide rational means for employing multiple rows or columns of 
the FRF matrix to improve modal vector estimates. 

(b) The use of synthesized frequency response functions rather 
than original experimental FRF's has two advantages: (1) far less memory 
is required to store the synthesis parameters than is required to store 
a complete FRF, and (2) the analysis bandwidth can be reduced to obtain 
better estimates of the modal vectors. 

(c) Much additional research is needed^ particularly in three areas: 

(1) application of modal tuning to testing of more complex structures, 

(2) use of synthesized FRF's for reference points where no actual shakers 
were located, and (3) use of residuals where modal density is very large. 

7. Appendix A - Single Column/Row Methods of Modal Parameter Identification 
References 2 through 4, and many other papers, describe the basis 

for using a single column or row of the FRF matrix to estimate modal 
parameters. It is assumed that the motion of the physical system can be 
adequately described by a set of simultaneous second-order linear differen- 
tial equations of the form 
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Mx(t) + Cx(t) + Kx(t) == f(t) (A-1) 

where 

x(t) = system displacement vector (nxl) 
f(t) « force vector 
M = mass matrix 
C = damping matrix 
K = stiffness matrix 

It is convenient to take the Laplace transform of Equation (A-1) giving 
B(s) X(s) = F(s) (A-2) 

where 

Laplace transform of displacement vector 
Laplace transform of force vector 
Ms^ + Cs + K 

(complex) Laplace variable 
Equation (A-2) can also be written in the form 

X(s) = H(s)F(s) (A-3) 

(A-4) 


where 

H(s) = (Bfs))"’ 


X(s) = 
F(s) = 
B(s) = 
s = 
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H(s) is referred to as the transfer matrix. H(s®io)) is referred to as 
the frequency response function (FRF) matrix. 

The element h^j(s) of H(s) can be written 


h,j(s) 


u e2n-2 . . ^2n-l . . . „ ^ . 

M ^2^ + ... + bg^^is -t- 

' det(B(s)) 


(A-5) 


For an n-dimensional system there will be 2n roots of the equation 


det(B(s)) = 0 


(A-6) 


and if these roots are distinct, then H(s) can always be written in 
the partial fraction form 

2n a,. 




(A-7) 


Where 


P|^ = kth root of Equation (A-6) 
aj^ = residue matrix for the kth root 

When the system is subcritically damped, as are practically all structures, 
the roots P|^ occur in complex conjugate pairs, which can thus be 
written in the form 

p« = -Ok ■*“k ■ Pk “ -“k - ’"k 


(A-8) 
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Modal vectors, U|^ , are defined as solutions to the homogeneous 
equation 


B(P|^)U|^-0 (A-9) 

and, in terms of these modal vectors, the transfer matrix can be written 
in the form 


or 


H(s) 


k=i =-i>k 



(A-IO) 


(A-11) 


From Equation (A-10) or Equation (A-11) it may be seen that each row and 
each column of the transfer matrix (and, hence, the FRF matrix) contains 
each modal vector multiplied by a component of itself. Thus, if the effect 
of a single root could be isolated from the effect of all the other roots, 
the mode U|^ should be identifiable from any row or column of the result- 
ing H|^(s) , where 

T 

H,(.) = (A-12) 

except those rows and columns corresponding to null elements of U|^ . 
References 2 through 4 describe various methods for attempting to isolate 
the contribution, H|^ , of a single mode in order to identify the elements 

k * 


of the modal vector u 
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Table 1. Effect of Subsystem properties on System Modes of a Weakly 
Coupled 2D0F System. 
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E-stimated Roots 

H^l 




Root 

Frequency 

Damping 

/Amplitude 

Phase 

1 

115.2 

0.4707E‘*01 

32.19 

0.1799E-05 

2 

115.0 

0.9298E-02 

2.799 

-0.1945E-05 

3 

119.0 

0.2640E-02 

104.4 

1.684 

4 

119.5 

0.5844E-02 

11.33 

0.2844 

5 

124.8 

0.2241E-01 

10.46 

-3.142 

6 

124.8 

0.5535E-02 

0.2597E-02 

0.7300E-04 


Estimated Roots 

*^21 




Root 

Frequency 

vamping 

Amplitude 

Phase 

1 

115.1 

0.3076E-01 

20.96 

-0.2731E-05 

2 

115.0 

0.6774E-02 

1.287 

0.2747E-05 

3 

119.0 

0.2638E-02 

105.5 

1.692 

4 

119.4 

0.6918E-02 

12.86 

0.7228E-01 

5 

124.8 

0J989E-01 

8.361 

3.142 

6 

124.8 

0.4I53E-02 

0.4751E-02 

-0.2091E-03 


Estimated Roots 

«51 




Root 

Frequency 

Damping 

Amp! itude 

Phase 

1 

115.0 

0.8869E-02 

0.9159E-01 

0.1025E-05 

2 

119.0 

0.2248E-02 

21.34 

1.479 

3 

119.6 

0.2688E-02 

21 .75 

-1 .693 

4 

121.8 

0.2591 E-01 

0.81 92 

-1 .075 

5 

124.8 

0.2374E-01 

0.9778 

-0.2076E-06 


Estimated Roots 

«61 




Root 

Frequency 

Dampi ng 

Amp! i tude 

Phase 

1 

115.0 

0.1177E-01 

0.1757 

0.1314E-04 

2 

119.1 

0.2258E-02 

22.18 

1 .438 

3 

119.1 

0.2723E-01 

1 .053 

0.2236 

4 

119.6 

0.2839E-02 

22.71 

-1 .742 

5 

124.8 

0.1433E-01 

0.1876 

-0.2673E-05 


Table 2. MDOF Cu»"ve Fit Parameters for FRF's of Dual Beam Model 
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Estimated Roots 





Root 

Frequency 

Damping 

Amplitude 

Phase 

1 

115.2 

0.4602E-01 

2.207 

-0.3383E-05 

2 

115.1 

0.3774E-01 

0.131 5E-06 

-2,032 

3 

115.0 

0.3210E-02 

0.4392E-02 

0.2486E-03 

4 

118.9 

0.2873E-02 

12.55 

1.589 

5 

120.0 

0.2840E-02 

13.00 

-1.568 

6 

124.8 

0.1538E-01 

0.2645 

-0.3002E-06 


Estimated Roots 

^26 




Root 

Frequency 

Damping 

Amplitude 

Phase 

1 

115.3 

0.6888E-01 

1 .838 

0.4750E-06 

2 

115.0 

0.3831E-02 

0.3873Er01 

-0.1847E-05 

3 

M 

118.9 

•1 OA A 

0.2962E-02 

A AAe7f- AA 

0.1951E-01 

12.19 

1 A te 

1.589 

•1 A 

t 

5 

1 CV.U 
124.8 

i C • / 9 

0.5839E-07 

** 1 » ^ f C». 

-0.8005E-01 

6 

124.8 

0.1038E-01 

0.1001 

-0.3857E-06 


Estimated Roots 

^56 




Root 

Frequency 

Dampi ng 

Amp! Itude 

Phase 

1 

115.0 

0.1351E-01 

5.757 

O.lOlOE-05 

2 

118.8 

0.3026E-02 

105.2 

1.591 

3 

119.7 

0.951 4E-02 

0.3043E-02 

-1 .827 

4 

120.5 

0.2598E-01 

17.09 

-0.9058 

5 

124.8 

0.7826E-02 

1.399 

-3.142 


Estimated Roots 

*^66 




Root 

Frequency 

Dampi ng 

Amp! Itude 

Phase 

1 

115.0 

0.9683E-02 

2.987 

0.3488E-05 

2 

118.8 

0.301 2E-02 

104.0 

1.592 

3 

119.4 

0.2926E-01 

22.91 

-0.7715 

4 

119.8 

0.3661 E-02 

0.7704E-01 

-2.299 

5 

124.8 

0.9513E-02 

2.130 

3.142 


Table 2.{Cont.) MDOF Curve Fit Parameters for FRF's of Dual Beam Model 











































BW 

(Hz) (Hz 


10 119.023 


4.0 119.024 



Mode 2 


Amp! . Phase 


1.000 0 . 0 ° 

1.001 0 . 0 ° 

0.200 -34,6° 

0.193 -17.3° 



119.023 


Mode 2 


Pt, Amp! . Phase 


1.000 

1.003 


0.048 -41.9° 

0.043 -40.0° 















































Figure 2. Frequency Response Function H-,, for 2D0F Systems, 
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Figure 3. Weakly- Coupled Systems Identified by Klosterman 





Figure 4. Bode Plot of Drive Point FRF of System with Closely-Spaced-Frequency Modes 







4.985 Hz 



*Frequenc1es based on cross- FRF quadrature peaks. 



Figure 7. Superposition of Frequency Responses of Closely-Spaced- Frequency f-fodes 
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Acquire single-shaker FRF's. 
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